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ABSTRACT: The recent LHC discovery of a Higgs-like boson at 126 GcV has important conse- 
quences for SUSY, pushing the spectrum of strong-interacting supersymmctric particles to high 
energies, very difficult to probe at the LHC. This gives extra motivation to study the direct pro- 
duction of electroweak particles, as charginos and neutralinos, which are presently very poorly 
constrained. The aim of this work is to improve the analysis of chargino-neutralino pair production 
at LHC, focusing on the kinematics of the processes. We propose a new method based on the 
study of the poles of a certain kinematical variable. This complements other approaches, giving 
new information about the spectrum and improving the signal-to-background ratio. We illustrate 
the method in particular SUSY models, and show that working with the LHC at 100/fb luminosity 
one would be able to distinguish the SUSY signal from the Standard Model background. 
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1. Introduction 

The recent discovery of the Higgs boson, with a mass around 126 GeV [1], [2], does not only 
have crucial importance by itself. It also has far-reaching consequences for well-motivated 
candidates for physics beyond the Standard Model (SM), such as Supersymmetry (SUSY), 
and in particular the Minimal Supersymmetric Standard Model (MSSM). 

As it is well-known, in the MSSM the tree-level Higgs-mass is bounded from above 
by the mass of the Z-boson, so large radiative corrections are needed in order to reconcile 
theory and experiment. These corrections go with (the logarithm of) the supersymmetric 
masses, in particular with the stop masses. Consequently, these masses must be rather 
high (well above 1 TeV unless there are unlikely conspiracies) in order to reproduce the 
experimental Higgs mass. This means that the MSSM is in trouble. First of all, such 
large supersymmetric masses are in tension with the naturalness of the electroweak break- 
ing. But, beyond naturalness considerations, large supersymmetric masses imply that the 
MSSM is going to be difficult to test at the LHC. 

In constrained versions of the MSSM, such as the popular constrained MSSM (CMSSM), 
the stop masses are directly related to the universal scalar mass, mo, and gaugino mass, 
m 1 / 2 , at the high scale. So large stop masses imply also large squark and gluino masses. 
Actually, the lower bounds for these physical masses are beyond the present (and most 
probably the future) LHC reach [3, 4, 5, 6, 7, 8, 9]. In other words, even if the MSSM is 
true, mh — 126 GeV implies that it will be difficult to discover SUSY at the LHC studying 
the production and decay of squarks and gluinos, which has been so far the golden channel 
of potential SUSY production. Going beyond the CMSSM does not help much, as the 
values of the physical stop masses are directly related to the stop and gluino initial masses 
(due to the coupled renormalization-group running). Still, one could have heavy stops and 
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light first and second generation squarks, though this possibility is rather strange from the 
theoretical point of view. On the other hand, the production of supersymmetric particles 
with no strong interactions may become one of the most promising possibilities to detect 
SUSY, especially in non-constrained MSSM models. 

The latter is an important motivation to explore chargino-neutralino production and 
detection at the LHC. Actually, relaxing the severe universality of the CMSSM (in partic- 
ular for the gaugino sector), one can have relatively light charginos and neutralinos, while 
still having squarks and gluinos very heavy, e.g. going into the focus-point region. On 
the other hand, charginos and neutralinos can produce not only hadronic, but also purely 
leptonic final states, with distinctive features which make them very attractive. 

Another motivation comes from dark matter (DM) constraints. The last data of 
XENON100 in combination with the Higgs mass have narrowed enormously the MSSM 
candidates for DM, see e.g. refs.[10, 7]. Probably the most satisfactory scenario that sur- 
vives occurs when the lightest neutralino is almost pure Higgsino with a mass around 1 
TeV. But if gaugino masses are not universal it is also possible to have a much lighter 
(mostly bino) neutralino which annihilates in the early universe through a combination of 
Z— boson and Higgs funnels, see [10]. 

In summary, a scenario where all the supersymmetric particles are too heavy, except 
charginos and neutralinos, is plausible and has phenomenological motivations. Therefore it 
would be worthy to improve the present techniques to analyze the production and detection 
of chargino/neutralino pairs at the LHC; and this is the main motivation of this paper. 

In most of the cases the chargino-neutralino pair created is xfx2> i- e - the lightest 
chargino and the next-to-lightest neutralino. Some of the diagrams of production and 
decay are shown in Fig. 1. The chargino and the neutralino can decay in several ways, 
always giving an LSP (x?) at the end of each cascade. 

The study of chargino-neutralino pair production has been performed previously in 
e.g. [11], [12], [13] through the analysis of leptonic final states. The aim of our work is to 
improve those analyses, by proposing new strategies which are complementary and more 
efficient in some cases . Our analysis is by construction independent of the diagrams 
through which the chargino and the neutralino decay, since it is entirely based on the 
properties of the final states. It can also be applied to final states including hadrons, such 
as bbl. 

The outline of the paper is as follows. In section 2 we describe the present status of the 
analyses of possible chargino-neutralino production and detection, as done by ATLAS and 
CMS groups, and motivate the convenience of improving this kind of analysis. In section 3 
we describe our strategy, by proposing a kinematical variable which allows to obtain direct 
information about the SUSY spectrum. It is also very useful to concentrate signal-events, 
thus improving the S/B ratio, and thus the LHC potential for discovery of new physics. 
Section 4 is devoted to illustrate the use of this variable in concrete SUSY models. We 
finally conclude in section 5. 

For a general review on kinematical techniques, see [14]. 
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Figure 1: a) and b) Typical production processes of a pair X1X2 m the LHC. c) and d) Chargino 
decay modes, e) and f) Neutralino decay modes. 

2. State of the art of the analysis 

As shown in Fig.l, the neutralino X2 can decay through a charged sfermion, a Higgs or a 
Z— boson, depending on the kinematical availability and the x% composition. Likewise, the 
chargino xf can decay, among other channels, through a charged sfermion or a W— boson. 
However, as long as the sleptons are heavier than x% an d xfi these decay channels become 
very suppressed and the decays through on-shell Higgs, Z and W dominate. 

Ref. [11] contains a study of purely leptonic final states of the type t + £~l !± , where I 
and £' may be identical leptons. The authors perform separate analyses of two cases: 1) 
the invariant mass of identical opposite-sign leptons (m^+i-) does not reproduce the Z- 
boson mass, Mz, and 2) it does. They assume that the X2 decays through a slepton (case 
1) or a Z— boson (case 2). In both cases, they use the E™ 1SS variable to compare the actual 
experimental data with the SM background, finding no significant excess of events so far, 
once all the uncertainties are taken into account. The negative result is then interpreted 
as contour bounds in the parameter space of e.g. concrete simplified models. 

A very complete analysis was presented later in ref. [12]. Again, the authors focus on 
3-lepton final states, but using either the variable E™ 1SS or rri£+£- in combination with Mt 
(the transverse mass built with the momentum of the unpaired lepton) and E™ ss - One 
option gives better sensitivity than the other, depending on the mass splittings between 
X± and its respective mothers. A second analysis was performed for 2-lepton final states, 
considering the possibility that one of the 3 leptons produced by X\ or X2 ma y be lost or 
does not pass the kinematical cuts. In a last analysis they considered that both xf and x\ 
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may decay through on-shell vector bosons, giving 21 + 2j in the final state, for which the 
SM background has not intrinsic E™ 1SS . Using these techniques they were able to provide 
contour bounds on X\ > X2 an d X? masses for particular simplified models. 

A quite efficient analysis of the decay through sleptons was presented in ref. [13] (based 
on [15] and [16]), where the authors used the mg+i~ variable. An obvious advantage of 
this choice us that rri£+£- is Lorentz invariant and thus the analysis is fully valid in the 
(boosted) LAB frame. It can be shown that, when the intermediate sleptons are produced 
on-shell, the histogram on rri£+£- has an edge for rri£+£- \ edge = /(M^o, M x , m^), where / is 
a certain combination of the X2>Xi and t masses. This potentially gives a very distinctive 
experimental signal. When the sleptons are produced off-shell, the authors use a strategy 
based on the study of end-points in m,£+£-. This case is much less efficient, since by 
definition one has poorer statistics and a more difficult-to-control background. 

In this paper we propose the use of a variable whose distinctive feature is to concentrate 
the signal events around a peak, not relying on end-points. It also provides direct infor- 
mation about the spectrum. Besides, it can be applied without any assumption or guess 
about the decay mode that takes place (through Z— boson, sleptons, Higgs or whatever). 
We describe the idea in the next section. 

3. Our Strategy 

3.1 The visible transverse energy 

In contrast with the previous analyses, our strategy is purely kinematical and based only 
on the characteristics of the initial (%°) and the final states so it is independent of 

the channel through which the X2 decays. As it is based just on the chain, the analysis 
can be easily combined with other analyses that use partial information from both the X2 
and the Xi chains. It can also be applied to other processes where one or two x\ states 
are produced. 

We initially work in the reference frame in which x% is at rest, which we call CM (do 
not confuse with the center-of mass of the partonic collision) and consider the St variable, 
defined as 

S T = E T + E T (3.1) 

where E T , Ej, are the transverse energies of the the visible system (e.g. v = £ + £~) and the 
missing system, 

E T = y/M* + (^) 2 , E T = ^Ml + (p£)2 (3.2) 

Here M v and M x are the invariant masses of the visible system and the LSP (X?) and hats 
denote CM quantities everywhere. Of course, (p T ) 2 = (p T ) 2 - 

As it has been discussed in ref. [17], the histogram of events displayed in the St 
variable has a pole at the mass of xi 

£t\ po1c = Ecu = M X2 (3.3) 

Besides, for St > M2 there are no events, so the histogram has an sharp edge at the pole. 
Hence, potentially, the St variable can give a very distinctive signal, well separated from 
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the background, providing in addition direct information about the SUSY spectrum. Note 
also that £t is invariant under longitudinal boosts. There are however some problems. 
First, even working at the CM frame, i.e. assuming that X2 was produced with vanishing 
transverse momentum, we cannot measure the invisible transverse energy, Ej,, due to the 
uncertainty on the value of M x . Actually, we have checked that in general M x ~ is not a 
good approximation. Second, the x% neutralino is usually produced with a non-vanishing 
transverse momentum. We postpone the second issue to the next subsection and focus now 
on the first one. 

In order to avoid the dependence on unknown quantities, a good strategy is to work 
just with the visible transverse energy, Ej,, defined in eq.(3.2). It can be easily checked 
that the pole in the £t variable, eq.(3.3), translates into a pole in Ej,, 



polo 2M- 



[M 2 X2 -M 2 x + Ml] (3.4) 



\2 



In general, the mass of the visible system, M v , can change from event to event. However, 
if X2 decays through a Higgs (very common case) or through a Z— boson, all the events are 
concentrated around M v = (Mz)- An example of this kind, which shows clearly the 
edge and pole in the CM frame can be seen in the plot a) of Fig. 2, to be discussed below in 
more detail. When X2 decays through sleptons, we can select a fraction of the events with 
similar M v , and check that in CM a pole appears as in (3.4). This decreases the statistics 
but not in a dramatic way. Actually, in cases where M 2 <C M^ 2 , one could do a histogram 
with all the events, since the pole appears then around Ej, ~ ^ — (M£ 2 ~~ ^x) m 
Next, we study what happens when we go from the CM to the actual LAB frame. 

3.2 From CM to LAB 

The clean and sharp behavior of the Ej, variable in CM becomes of course less keen-edged 
when passing to the LAB frame. This is illustrated by the plots of Fig. 2. They correspond 
to a CMSSM model with m = 500 GeV, M 1/2 = 700 GeV, tan /3 = 10, A = 0, and /j, > 0. 
The associated relevant spectrum for us is M X2 ~ 554 GeV and M x ~ 294 GeV. The X2 
neutralino decays almost entirely through an on-shell Higgs. In the CM frame the pole for 
Ej,, eq. (3.4), is at ~ 210 GeV. We have simulated the events using Pythia, with a center- 
of-mass energy y/s = 7 TeV and a luminosity of 20/fb. The upper (lower) plots of Fig. 2 
correspond to events where the Higgs decayed into rf (a pair of b-jets, jbjb)- The left (right) 
plots are obtained in the CM (LAB) frame. We emphasize that this example is just for 
illustrative purposes. For this reason we have not incorporated a realistic r-reconstruction, 
as well as b-jet identification. 

As expected, the CM histograms show a very clear and sharp peak (pole) around 
200-220 GeV, followed by an edge. In the case of the jbjb histogram there appear events 
beyond the edge due to the limited efficiency in the jet reconstruction. The corresponding 
LAB histograms show two (non-dramatic) differences. First, a logical spreading around 
the maximum; and second, a slight shift of the peak towards smaller values of Ej-. This 
shift is almost invisible for the jbjb histogram. All these effects can be explained, estimated 
in a semi-analytical way and kept under control, as we discuss below. 
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Figure 2: Histogram (number of events) of in the cases where v is a) a rf pair in the CM 
frame; b) a rf pair in the LAB frame; c) a jbjb pair in the CM frame, and d) a jbjb pair m the LAB 
frame. The details of the SUSY point correspondent to this histograms are described in the text. 



We recall that the transverse variables we are using are only affected by transverse 
boosts when passing from CM to LAB. Longitudinal boosts are irrelevant for this analysis. 
There are two sources of transverse boosts. First, the partonic collision will not occur 
in general in its center-of-mass, i.e. it will have a non-vanishing net transverse momen- 
tum, typically due to initial state radiation. Second (and more importantly) even at the 
center-of-mass of the partonic collision, can be produced with non-vanishing transverse 
momentum (and opposite to that of xf) ■ Next we estimate the change in the EV^ of the 
events at the pole due to the non- vanishing transverse momentum of \\ in the LAB. 

We will distinguish two perpendicular directions in the transverse plane: the direction 
along p%? (||), and the perpendicular to it (_L). Due to the boost in the || direction (we will 
ignore the boost in the longitudinal direction, which is irrelevant here), the visible energy 
changes (from CM to LAB) as 

E v ^ E v = jE v — (3 j (pt)\\ (3.5) 

where (Px)\\ is the component of the visible 3- momentum in the || direction. As usual, 7 
and f3 are the parameters of the Lorentz transformation, satisfying 

^ Jj P- < 3 - 6 > 

1VI X2 

Now, for the events at the pole E^ = E v in the CM frame [17]. This holds after the 
transverse boost. Hence, when going to LAB, for those events the transverse energy of 
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those events changes as in eq.(3.5): 



which represents a shift 



£*->££= 7 £*-0 7 (fcy),l (3.7) 



A££ = (7- 1)^-07 (£r)|| • (3-8) 



Let us estimate the size of this shift. Since for a non-relativistic boost, 7 — 1 ~ ^/3 2 7 2 , 
the second term in the r.h.s. of (3.8) will normally be the dominant one, because it scales 
with f3 instead of /3 2 . However this term is sometimes positive and sometimes negative 
(depending on the sign of {$), whereas the first term, which scales with /3 2 , is always 
positive. On the other hand, in average, (py)|| ~ ~^Pt ■ Using eq.(3.4) we can express 
as a combination of the masses involved in the system, 

(^) 2 | polc = Tj^- [< - Ml + M 2 ] 2 - Ml . (3.9) 

This gets simplified when M 2 <C M^ 2 , e.g. when the decay occurs via Higgs or Z— boson 
and also in the other cases if we select events with relatively small M v . Therefore, at the 
pole 



pole 



which substituted back in eq.(3.8) gives 



AE* T ~ Q/3 2 7 2 - -L/3 7 ) E% . (3.11) 

We recall that this is expression is valid for events lying at the pole, which of course are 
especially interesting as they provide the maximum in the histogram of the signal. Now 
for a numerical evaluation of AE^ we need to estimate /3 7 . It can be empirically checked 
that typically the transverse momentum of \% * s i n average 

W\-\M X2 (3.12) 

and this value gets slightly smaller as M X2 increases. Actually this is consistent with 
expectations. First, the production of X1X2 through a W in 5-channel (see Fig. la) is 
dominated by the resonance. Even though the W will be normally off-shell, the penalisation 
for large momenta of X\X% ls important, although this reason becomes weaker for M\y <C 
M X2 + M^±. Second, and more importantly, the PDFs penalize large energies. Since 
the final state particles are colourless (in this case charginos and neutralinos) with a net 
electric charge, they are essentially produced via quarks/antiquarks (and not gluons). This 
is significantly affected by the PDFs of the anti-quarks, which tend to prefer lower momenta 
rather than higher ones. As a consequence of these two effects, the energies of the particles 
at the final state are expected to be very close to their mass. This is consistent with the 
relation (3.12), which implies that the energy of the particle produced differs less than 5% 
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Figure 3: Scatter plot of different SUSY models showing the correlation between E^\ po \ e and the 
peak value of the variable E^, reconstructed event by event for each point. This is done for the 
cases: a) rf pair in the CM frame; b) rf pair in the LAB frame; c) jbjb pah in the CM frame, and 
d) jbjb pair in the LAB frame. For the plot correspondent to taus, the green points show SUSY 
models where the neutralino \2 decayed through f instead of H. The straight line is the theoretical 
prediction in the CM frame. 



from its mass. In consequence, from eq.(3.6, 3.12) typically /?7 ~ 1/3. Then eq. (3.11) 
becomes 

where we have already incorporated the fact that f3 has a random sign. This implies a shift 
in the range [— 0.18-Ey, 0.29E^] in the visible transverse energy of the pole events when 
going from CM to LAB. Finally, we have to evaluate how this modifies the position of the 
peak. Notice that the pole events shifted positively (to the right of the histogram) , will not 
produce any new global maximum as they fall in a region of EVp where there were almost no 
events. On the other hand the pole events shifted negatively (to the left of the histogram) 
will populate a region where there were already events. In addition, that region gets also 
populated by positive shifts of non-pole events, with lower Ej, values. Consequently, we 
expect a maximum in a value of EJ^ which is approximately 18% smaller than the value of 
the pole in CM, see eq.(3.4). 

Fig. 3 shows scatter plots from a large set of SUSY models, showing the correla- 
tion between the theoretical £^| po i e , given by eq.(3.4), and the actual peak value of the 
EVp— histogram, reconstructed for each SUSY model. The upper (lower) plots correspond 
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to histograms of events where the X2 neutralinos decayed into a rf- (jbjb-) pairs. We 
remark again that these plots are shown for illustrative purposes only, so no attempt of 
realistic identification of the r and b— jets states is done at this level 2 . The left (right) 
plots correspond to histograms in the CM (LAB) system. The CM plots show the perfect 
correlation of the histogram maximum and the theoretical pole (3.4). The LAB plots show 
the two effects expected: a certain spreading and a net shift of the average maximum with 
respect to the CM prediction (straight lines). More precisely, the results for rf events are 
consistent with the previously discussed AE^ ~ — 18%E^, expectation. It is remarkable, 
however, that for jbjb events, the LAB results are more symmetrically distributed around 
the CM prediction. This can be easily understood by recalling that for jbjb events the 
CM histogram does not have an end-point at the pole because of the limited efficiency in 
the jet reconstruction, see plot c) of Fig. 2. In consequence the discussion after eq.(3.11) 
gets now modified since the pole jbjb~ events shifted positively fall in a region of Ej, where 
there were already events. In consequence, the limited efficiency in the jet reconstruction 
funnily makes the position of the maximum more stable than for the leptonic histogram; 
thus the symmetric spreading of the LAB scatter-plot around the CM prediction. 

One of the most interesting features of the Ej, variable is that it concentrates the signal 
events around a maximum (determined by the supersymmetric spectrum), which does not 
happen for the background events (these get concentrated around Mz, see below). This 
behaviour is very helpful to separate the signal events from the background events, thus 
improving the signal/background (S/B) ratio in the search of new physics at LHC (even 
before using this variable to extract information about the SUSY spectrum). It should 
be noticed here that the background events (typically SM production of W Z) behave as 
if the "neutralino" X2 na, d exactly the mass of the Z-boson and the "neutralino" Xi na d 
zero-mass. Thus for the background events the peak in the Ej, variable lies at ~ Mz- 

For some SUSY models, however, the position of the maximum of the signal -determined 
by eq.(3.4)- can be close to Mz, i.e. the peak of E^ for background events. Then, the 
strategy to detect the existence of new-physics events can be further improved by plotting 
the p?r variable rather than E^ (of course, both variables are equivalent since they are 
unambiguously related through E^, = \JM^ + (p^) 2 , where M v is a measurable quantity). 
Note that for the background events, the peak in pj, is at zero (in practice, not far from 
the lower kinematical cut used), while for the signal events is given by eq.(3.9), and lies 
normally at some non-trivial value (even if the corresponding Ej, is close to Mz). In the 
next section we will see some explicit examples where the use of the pVp variable is very 
convenient to show the existence of new physics in the first place. Once the new physics 
is detected, the value of pj* at the maximum of the signal can be related to the SUSY 
spectrum through eq.(3.9) or equivalently through eq.(3.4). 

Next we explore these features in further detail. 

4. Testing the efficiency of and p v T in concrete SUSY models 

As we have described before, the chargino-neutralino pair production can be studied 
2 In this sense, a more realistic analysis is performed in the next subsection 
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through the use of the visible transverse energy, E T , regardless of the way the neutralino 
X 2 has decayed. So here we only care about initial and final states. 

We will simulate LHC signals with y/s =14 TeV and luminosity of 100/fb, using the 
codes SuSpect , HDECAY and SDECAY 3 as well as S0FTSUSY[22] for the spectrum calculators, 
and MadGraph /MadEvent [23] and PYTHIA[24] for the event simulation. We focus on events 
with 3 leptons + missing transverse momentum, and apply the following general cuts 

• the existence of at least two identical, opposite signed leptons 

• the 3rd hardest lepton having Et > 10 GeV 

• At least an electron (muon) with Et >25 GeV (pr >20 GeV) 

• 3 leptons with pt > 20 GeV and 77 < 2.47 for electrons (77 <2.4 for muons) 

• E T liss > 50 GeV 

• The transverse mass of the 2nd chain Mt > 90 GeV, 4 

• Jets p T < 20 GeV 

• PGS ATLAS cuts [23] . 

The jet reconstruction is performed using the anti-fcr algorithm with AR = 0.4. 

We illustrate our strategy by working in the context of the MSSM. Specifically we choose 

the following SUSY models (defined at the Mz scale) 



Model 1 
Model 2 
Model 3 



Mi ~ 99 GeV, M 2 ~ 183 GeV, /i ~ 705 GeV, tan/3 = 10 , 
Mi ~ 47 GeV, M 2 ~ 244 GeV, /J ~ -515 GeV, tan /3 = 19 , 
Mi ~ 93 GeV, M 2 ~ 405 GeV, n = -5372 GeV, tan p = 50 



(4.1) 



where Mi , M 2 are the bino and wino mass parameters at low-energy, and tan (3 = (H u ) /(H^) 
is the ratio between the VEVs of the two Higgs doublets. We do not specify the values of 
gluino and sfermion masses, which are assumed to be heavy (note that we are not imposing 
gaugino-mass unification) . In the three models we study the signal of x 2 decaying through 
a Z— boson, which is the dominant one for all of them (more details below). The dominant 
background is WZ and VFZ+jet production in the SM. 

The strategy of the analysis for the three models is the following: 

• We reconstruct the possible invariant masses M v of a pair of identical leptons, select- 
ing the events where a pair of opposite-sign leptons has an invariant mass close to 
M z 

M V = M Z ± 10 GeV . 

3 all included in a single package called SUSY-HIT [18, 19, 20, 21]. 
4 Computed with the unpaired lepton, see discussion below. 
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This pair is identified as a daughter of a Z, thus the third lepton, £', should come 
from a W. This holds for both the signal and the SM background. 

It is important to note that most of the WZ background events are actually removed 
thanks to the Mt > 90 GeV requirement in the above list of cuts. The reason 
is the following. Taking into account that for the events of the background the 
missing momentum can be identified with the neutrino, one constructs the associated 
transverse mass, which satisfies the following inequality 

(M^f = (£# + E" T f - (p/ + p T n 2 < Myy , (4.2) 

where £' is the charged lepton coming from the W (i.e. the unpaired lepton). Hence, 
identifying p™ lss with p^ and discarding events with Mj, v < M\y would in principle 
remove the whole WZ background. Although, due to the inefficiency in the recon- 
struction of the missing piece, some background events may violate in practice the 
inequality (4.2). This turns out to be an extremely efficient constraint to improve 
the S/B ratio. 

• With the surviving events, we reconstruct the Ej,, p\ variables associated to the 
lepton pair daughter of the Z. 

• We simulate the background taking into account only the contribution from WZ and 
WZ+jet production, which is by far the dominant one in the regions of interest (see 
e.g. [11]). We have checked that the results of our simulations are in agreement with 
[11]. 

• Finally, we construct the Ej,, pVp histograms to show the signal over the background. 
Next we expound the results obtained for the three models considered, in a separate 

way. 

Model 1 

From the initial parameters listed in eq.(4.1), we obtain the relevant supersymmetric 
spectrum for the analysis, which reads 

M + ~ 207 GeV, M q ~ 203 GeV, M Y o ~ 107 GeV. (4.3) 

A,i 

Since the staus are heavy, a priori the preferred decay channel of \% i s through the lightest 
Higgs, rrih = 126 GeV. However, since M X2 — M x < mh, this channel is suppressed and the 
favourite decay-channel turns out to be through an on-shell Z. As mentioned above, the 
main background for this final-state topology is the SM production of a W/Z pair. 

Fig. 4 shows the histogram corresponding to the variable Ej,, considering the back- 
ground plus the signal with a luminosity of 100/fb, which will be reached in the near future. 
The SUSY signal is clearly visible. The position of its peak can be used to extract infor- 
mation about the SUSY spectrum in the way discussed in sect. 3. More precisely, for this 
model the prediction of eq.(3.4) for the visible transverse energy at the pole of the signal 
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Figure 4: Histogram (in number of events) of Ej, for the parameters corresponding to Model 1 
(see text), taking into account the signal plus the dominant background (W/Z production). 




ooooooooooooooo 
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Figure 5: The same as in fig. 4, but now the histogram is on the variable pVp for: left) background 
+ signal, right) only signal. 



is Ej, ~ 94 GeV, which is consistent with the bin of the histogram corresponding to 

pole 

the maximum signal, i.e. 90 GeV — 110 GeV. 

This case illustrates the possibility that the peak of the EVp— histogram for the signal 
lies close to the peak for the background, i.e. Mz- Note from Fig. 4 that both peaks are 
at the same bin, as expected. 

As discussed at the end of sect. 3, this shortcoming can be improved by plotting p?p 
instead of Ej,. The result is shown in Fig.5(left). 

Although the improving is not dramatic, it is clearly visible. Note that the peak of the 
signal (the bin centered at 50 GeV) is now displaced with respect to the background one 
(the bin at 30 GeV). Fig.5(right) shows the Pj^— histogram of just the signal events. Of 
course, it is not something one can realistically obtain from the experiment but we present 
it in order to show the shape of the signal events, gathered around the (theoretically 
predicted) maximum. 
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Model 2 

In this case the relevant supersymmetric spectrum reads 

M + ~ 273 GeV, M o ~ 273 GeV, M v o ~ 53 GeV. (4.4) 

Xi A2 Al y ' 




oooooccooooo 



E\ (GeV) 

Figure 6: Histogram (in number of events) of Ej, for the parameters corresponding to Model 2 
(see text), taking into account the signal plus the dominant background (W/Z production). 
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Figure 7: The same as in fig. 6, but now the histogram is on the variable for: left) background 
+ signal, right) only signal. 

X2 an d Xi are almost pure bino and wino respectively. Now the decay of X2 through 
an on-shell Higgs is kinematically allowed, so one would expect it to be the preferred decay 
channel. However, the branching ratio of the decays through a Higgs or through a Z— boson 
are 39% and 61% respectively. The reason is the following. Due to the relatively high value 
of tan/3, the x\ an d X? neutralinos are essentially gauginos with a small HQ component 
and a very small H® component. On the other hand, the physical Higgs-boson is essentially 
H® for the same reason. Then the X2 decay through a Higgs occurs thanks to those very 
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small H® components and gets suppressed, while the decay through a Z— boson may occur 
also through the not-so-small H® components of the neutralinos. 

Fig. 6 shows the histogram corresponding to the variable E^, containing the back- 
ground plus the signal, again with a luminosity of 100/fb. The SUSY signal is of course 
weaker than in Model 1, since the supersymmetric masses are larger. On the other hand, it 
is nice that now the peak of the signal does not coincide with the peak of the background. 
In this case eq.(3.4) gives Ej, ~ 146 GeV, consistent with what the signal peak visible 

pole 

at Fig. 6. 

As for Model 1, in this case the separation between the signal and background peaks 
is more efficiently achieved using the variable pVp rather than Ej,, as it is shown in Fig. 7. 

Model 3 

Finally, for Model 3 the spectrum reads 



M + ~ 405.3 GeV, M v o 

Xl A2 



405.3 GeV, M v o 



94 GeV. 



(4.5) 



Again here %2 an d X? are almost purely wino and bino respectively, with very small 
Higgsino components of C(< 10 -2 ). Also in this case the neutralinos contain much more 
H® than H®, leading to branching ratios BR(^2 —> xV 1 ) and BR(x° — > Xi^) to be 4% and 
96% respectively, for similar reasons as before. Since now tan /3 is much larger the Higgs 
channel is even more suppressed. 

Fig. 8 and Fig. 9 show the histograms in the E^ and the p^ variables respectively. 
They show similar features as for Model 2. The difference of course is that the signal is 
now quite small due the large supersymmetric masses. Still it is visible, especially for the 
Pj, variable. Since the S/B ratio is quite good around the p^-peak of the signal, the later 
would be clearly visible with larger luminosities (say 300/fb). 



£ = 100/fb 



LOO 




E V T (GeV) 



Figure 8: Histogram (in number of events) of Ej, for the parameters corresponding to Model 3 
(see text), taking into account the signal plus the dominant background (W/Z production). 
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L = 100/fb £ = 100/fb 
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Figure 9: The same as in fig. 8, but now the histogram is on the variable pVp for: left) background 
+ signal, right) only signal. 

5. Conclusions 

The study of chargino-neutralino pair production at the LHC finds nowadays an unprece- 
dented motivation, since the detection of those particles could well represent the best 
(maybe unique) hope of discovering supersymmetry at the LHC. Different strategies have 
been used so far by the experimental groups ATLAS and CMS to explore this process. They 
have been mainly focused on the analysis of 3-leptons plus missing-energy final states using 
standard kinematical variables, such as jrp and the invariant mass rri£+£-. 

In this work, we have presented a new -purely kinematical- method based on the 
variable Ej, (visible transverse energy) , or equivalently (visible transverse momentum) , 
which presents some very useful features. First of all, the histogram in Ej, (in the frame 
where the decaying neutralino is at rest) has a pole. This translates into a peak in the 
actual experimental distribution. Consequently, this is a very robust feature against poor 
statistics. In addition, this concentration of signal-events around the maximum does not 
occur for the background-events, thus improving the S/B ratio and the discovery potential 
of new physics. Finally, the peak is correlated in a well-defined way to a precise combination 
of SUSY masses, namely 

£tU = ^ - M l + M "] (5-1) 

where M X2 , M Xl and M v are the masses of the decaying neutralino \2-> the lightest neu- 
tralino (LSP) and the visible system to which decays to (typically two leptons or two 
b-jets). 

Of course, when passing from the CM to the LAB system the maximum becomes less 
sharp and the position of the maximum is slightly shifted. However, these effects are not 
dramatic and, besides, can be well estimated in a semi-analytical way and kept under 
control. 

We have illustrated these facts performing a realistic analysis of particular SUSY mod- 
els, where the \2 neutralino decays mainly through a Z boson. We have focussed in events 
where the final state consists of 3 leptons plus missing energy, showing that the SUSY 
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signal in the Ej, variable could be very well detected at the LHC running with 14 TeV and 
a luminosity of 100/fb, which is attainable in the near future. 

Let us finally stress that the same method can be applied to any possible decay channel 
of X2- through squark, sleptons or a Higgs. The latter case is actually the most frequent 
one for MSSM models. Then, the most probable final state contains 2 b-jets, 1 lepton (or 
jet) and missing energy. The h — > bb channel has the additional advantage of having large 
branching ratios compared to the Z — > 11 case. On the other hand, the topology jbjb^ is 
much less clean and difficult to reconstruct properly. Thus it requires a separate study 
which is out of the scope of the present paper and will be the subject of a future research 
work. 
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